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Abstract 

The boson-conserving one-nucleon transfer operator in the interacting boson 
model (IBA) is reanalyzed. Extra terms are added to the usual form used for 
that operator. These new terms change generalized seniority by one unit, as 
the ones considered up to now. The results obtained using the new form for 
the transfer operator are compared with those obtained with the traditional 
form in a simple case involving the pseudo-spin Bose- Fermi symmetry U^{6)0 
U^{12) in its U^^{5) U^{2) limit. Sizeable differences are found. These 
results are of relevance in the study of transfer reactions to check nuclear 
supersymmetry and in the description of /3-decay within IBA. 
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The interacting boson model (IBA) has been extensively used in order to describe spec- 
troscopic properties of even-even and odd-even nuclei One of the observables studied 
within IBA is the stripping and pick-up single-nucleon transfer spectroscopic intensity be- 
tween even-even and odd-even nuclei. The form of the one-nucleon transfer operator in IBA 
is an open problem of current interest for several reasons: i) some calculations indicate that 
the form commonly accepted is insufficient for detailed calculations 0], ii) it is one of the 
key points for checking nuclear supersymmetry IQ,^, and iii) it is a basic ingredient for the 
study of /3-decay since in IBA it is described as a composed process of a neutron (proton) 
pick-up followed by a proton (neutron) stripping Several forms for the one-nucleon 
transfer operator in IBA have been proposed in the past based on different microscopic 
approaches. The original form, proposed by Scholten and Dieperink [^|lOl! was obtained in 
the generalized seniority (v) scheme [|1T| by using the Otsuka-Arima-Iachello (OAI) mapping 



method |]T2| and the number approximation [|T3l. This is still the operator implemented in 



the standard IBA codes and the most commonly used. Alternatively, other forms based 
on the nuclear field theory techniques the use of RPA for quadrupole phonons and 
mapping the fermion space onto the fermion-quadrupole space |TH], and the generalized 



Holstein-Primakoff mapping method [0,|T3] have been proposed but there have been few 
actual calculations using these operators. 

In this paper we reanalyze the standard IBA boson-conserving one-nucleon transfer oper- 
ator proposed originally by Scholten and Dieperink and improve that operator by including 
extra terms which change seniority by one unit and are the next order in d — boson creation 
and annihilation operator. We show that this new type of terms give rise to sizeable changes 
in the one-nucleon transfer spectroscopic intensities and consequently conclude that they 
must be included in the standard IBA operator. 

The starting point is the IBA boson-conserving part of the one-nucleon transfer oper- 
ator derived semi-microscopically by Scholten and Dieperink P,p!0| using the generalized 
seniority formalism |TT]]. Consequently, it is only strictly valid for spherical nuclei where 
seniority is approximately a good quantum number. However, this operator has also been 
used with reasonable success for transfer reactions involving deformed nuclei Scholten 
and Dieperink constructed the Interacting Boson-Fermion Approximation (IBFA) image of 
the shell model single-nucleon creation operator, taking into account terms which change 
generalized seniority by one unit and using fermion states with v < 2. Under these approx- 
imations the transfer operator which does not change the number of bosons is 



where is the odd-nucleon creation operator in the IBFA boson-fermion space, Vj is the 
occupation probability of the single particle orbit j and Uj = ^1 — Vj. N is the number of 
bosons. The summation on j' runs over the valence single particle angular momenta, and 
the normalized quantities Pj'j are PJTO(| 

Pj'j = Qj'j (uj'Vj + Vjmj)/Np , (2) 

where 

Qfj = UW^'^Wj) (3) 
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and Nf3 is a normalization constant obtained from the condition 



T.iPn'Y = 1 (4) 



where j and j' are the valence single particle angular momenta. The operators dm , which 
behave appropriately under rotations, are dm = (— l)™'(i-m • 

As mentioned before, the expression (|l|) was obtained under the restrictions of changing 
generalized seniority by one unit and using fermion states with v < 2. The term [((i^ x 
dY'^'' X a],]^-* changes generalized seniority in one unit as [(s^ x d)^"^^ x aj,]^\ and does not 
appear in the traditional expression of the one-particle transfer operator. The reason is 
that the derivation of the corresponding coefficient implies the use of fermion states with 
V > 2. However, this kind of terms are obtained naturally when using nuclear field theory 
techniques |]T^, or the generalized Holstein-Primakoff mapping method [|T^. For this reason 



it is natural to wonder which is the infiuence of this kind of terms. Taking into account 
these terms, the new transfer operator reads Q 

^■i; = + E <i^jfH X d^y'^ X ^i^- (5) 

This is the operator to be used in IBFA if no distinction between protons and neutrons is 
made. When working with versions of the IBA-IBFA in which the proton-neutron degrees of 
freedom are treated explicitly, all the operators in the preceding equation will have an index 
TT (u) if the transferred nucleon is a proton (neutron). In this work the coefficients (pjj, will 
be obtained by requiring that the matrix elements of the fermion operator C] and the IBFA 



operator T^' between states of f = 2 and v = 3 in the corresponding spaces are equal. For 



that purpose, first one has to construct the v = 3 shell model states which, in general, will 
not be orthonormal. From this set of non-orthonormal fermion states a mapping procedure 
onto the v = 3 orthonormal boson states is needed. In this paper we use the so called 



democratic mapping ||18l. This mapping is based on the diagonalization of the shell model 
overlap matrix. In our case, the method starts from a set of n fermion states which are 
linearly independent but not necessarily orthonormal \F, (i); JM) {i = 1, . . . ,n) and a set 
of n orthonormal boson states \B, (i); JM) {i = 1, . . . ,n). The democratic mapping follows 
the following steps: 

1. First, one has to construct the initial fermion states \F, (i); JM) {i = 1, . . . ,n) and 
the orthonormal boson states \B, (i); JM) {i = 1, . . . ,n). 



We prefer to write the operator in this way but with the appropriate angular momenta recoupling 
it can be written as 
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2. Second, one has to construct the matrix G"^ with the overlaps of the fermionic states 
for each total angular momentum J, 

Qi = {F, it); JM\F, (k); JM) ; i, k = 1, . . . ,n. (6) 

3. The diagonalization of the matrix G'^ provides its eigenvectors and eigenvalues. Let 
us call C"^ the n x n square matrix containing the eigenvectors (in columns) fixed by 
the condition J^kiCfk)"^ = 1 and A"' the diagonal matrix containing the eigenvalues, 
A/. The following relations hold in matrix notation 

QJ .C^ = C-^ . (7) 

and 

{C^Y={CY- (8) 

The orthonormal eigenvectors are, 

1 " 

\F,k■,JM)^ = ^Y.Cf,k\F,{^)■,JM)■, k = l,...,n. (9) 



1=1 



4. The boson states are transformed to a new orthonormal set according to the same 
transformation as the fermion states 

n 

\B,k;JM)^ = Y.Cl,\B,{t);JM) ; k = l,...,n. (10) 

i=l 

It is easy to see that these new boson states are orthonormal using Eq. (|^). 

5. Finally, the mapping is done between the new fermion, Eq. (P), and boson, Eq. (p!OD, 
states 

\F, k; JM)i_ ^ \B, k; JM)^ ; = 1, . . . , n. (11) 



In the following, we proceed with these steps and derive semi-microscopically, in the 
same spirit as Scholten and Dieperink, the coefficients 0^^/. 

First, it is necessary to construct properly the '5 = 3 states in the shell model space. In 
order to get this set, we start with the following unnormalised shell model states 

|5^-^(Dj); JM) = [C] X D% (St)^-^IO), (12) 

which, in general, contain v = 1 and v = 3 components. Cj represents the creation operator 
of one nucleon in the single particle orbit j and the fermion pair creation operators S'^ and 
are [H, 
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s'=j:^J^ic]cjr (13) 

and 

^^ = E/5;'/|(C^jC7t)(^) (14) 

33' * 

where the structure coefficients aj and fi'^^i are, 

«. = \/¥^.' /5;.' = ^, (15) 
\ N UjUj' 

where Pjj' are defined in Eq. (j^), Qj = (2j + l)/2 and Qe = (2j + l)/2 is the effective 
shell degeneracy 0. 

Unnormalized v = 3 states are constructed subtracting the v = 1 components in ([T2|) in 
the following way 

3; S'^-'m; JM) = '^7^^^^'^' " = 1' JM), (16) 

/ (Af-i)!r(n,+i) 



where the denominator in the first term on the right-hand-side is introduced for later conve- 
nience in order to simplify the expressions of some quantities given below. \v = 1,N; JM) 
are properly normalized seniority v = 1 states. Xj is obtained requiring that the overlap 
between the states ([16|) and the states of seniority one is zero. The expression obtained is 



if J is equal to one of the allowed angular momenta, and zero otherwise. In the above 
derivation the number approximation [jl3[ (0|S'^(S'"'')^|0) = j^^^^j'^j^^^ has been used. 

Once we have the non-orthonormal fermion states, Eqs. ( [l6l) and (0), the matrix of the 
overlaps, Q'^ , can be constructed, obtaining 

Qi, = {€ = 3; S^'-'iDjrny, JM \v = 3; S^^^'iDjr); JM) 

= + 10 (f: 5,, E P,,,A,3^ { t 2 J } • (1^) 

Note that this expression tells us that the states (|l^) are already orthonormal when J does 
not coincide with one of the valence single particle angular momenta. For the general case 
in which there are several single particle angular momenta j that contribute to a given J one 
has to follow the procedure sketched above. First, the matrix B'^ is diagonalized obtaining its 
eigenvectors, C'^, and eigenvalues, A"^. Then, orthonormalized fermion and boson seniority 
three states are constructed as given in Eqs. @) and ([lOD 

1 " 

\F,v = 3,k-JM)^ = ^Y.CUv = 3-S''-\Dj,)- JM) ; k = 1, . . . ,n. (19) 

Aj, i=l 
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and 



\B,v = ?>,k-JM)^ = Y.Cl^\s''-\d3i); JM) ■ fc=l,...,n. (20) 



We are now in a position to map the f = 3 shell model states into the f = 3 boson-fermion 
states of the IBFA 

|F, v = ?,,k] JM)^ ^\B,v = 3, k; JM)±, k = l,...,n. (21) 

Since now the v = 3 states are known both in the fermion space and in the boson-fermion 
space, the coefficients 0/ will be obtained by requiring that the matrix elements of the 

fermion operator Cj^ and the IBFA operator Tj^ between states of v = 2 and -u = 3 in the 
corresponding spaces are equal, 

±{F,v = 3,i;J\\Cl\\S^'^D-2) =^ {B,v = 3, i; J\\T^^\\s^-^d;2) . (22) 

With this requirement the following expression for the coefficients 0^ ■ , is obtained 



Jk 



J 



^^u.-{c'-^-iCr)J, (23) 



where j = a/2j~+T and v A-^ is a diagonal matrix whose elements are the square root of 
of the corresponding ones in A'^ (the square root of the eigenvalues of G'^). The transfer 
operator given in Eq. (^ with the coefficients in Eq. ( [23| ) is the form proposed in this work 
for a particle-like transfer. For the case of a hole-like transferred nucleon the u's and v's 
have to be interchanged (in addition it should be noted that the reduced matrix elements 
given in Eq. (^) change sign if the transferred nucleon is hole-like [0). This new form for 



the boson-conserving one nucleon transfer operator has been implemented in the IBA codes 
| 20| which are available under request. 

We have performed calculations both in IBA-1 and IBA-2 with these new terms in order 
to show their relative importance and as a result sizeable differences are found. Below we 
present a simple example. 

The observable we calculate is the spectroscopic intensity defined as, 

/.(^.-^/) = l(^/ll^-^'ll^.)r- (24) 

Jj is a state in the even-even nucleus and Jf a state in the odd-even. It is easy to show that 
Ij{Ji — >■ Jf) = Ij{Jf — Ji) and this fact allows us not to specify which is the initial nucleus 
in the table presented below. 

As a simple example we present here one-nucleon particle-like transfer intensities between 
two nuclei which follow particular dynamical symmetries of the IBA-1 and IBFA-1. The 
even-even nucleus has just one boson, is the boson core for the odd-even nucleus. The odd- 
even nucleus has one boson plus one particle and is described within the pseudo-spin Bose- 
Fermi dynamical symmetry U^^{5)^U^{2) obtained from the product f/-^(6)®f/^(12), 
where the superscripts B and F stand for boson and fermion. In the odd-even nucleus, the 
odd particle has been allowed to occupy the j = 1/2,3/2 , and 5/2 single particle states. 
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Under these conditions, once the number of bosons is fixed, the wave functions are fixed by 
the symmetry, independently of the parameters used in the Hamiltonian. The occupation 
probabihties that appear in the transfer operator are not provided by the symmetry. In the 
example presented here the values used are, = 0-5, "^^3/2 = 0-25 and "^5^2 = which fulfil 

The states in the even-even core nucleus are identified by the labels of the irreducible 
representations of the group chain [/^(6) D U^{b) D O^iji) D 0-^(3): \[N]{nd){v)L) . In our 
case, with just one boson, there are only two states, the ground state |0b) =|[1](0)(0)0) with 
one s — boson and the excited state \2b) =| [1] (1)(1)2) with one d — boson. The states in the 
odd-even nucleus are labelled by the labels of the irreducible representations of the group 
chain f/^^(6) D f/^^(5) D 0^^(5) D 0^^(3) D Spin^^{?>): \[Nr,N2]{ni,n2){vuV2)L- J). 

Since the ground state of the even-even nucleus has no d — bosons the only term of the 
transfer operator that contributes to the transfers from this state to any state in the final 
odd-even nucleus is the a\ term. We have checked that our results in this case coincide with 

jk 



those obtained analytically in Ref. [^. The calculated spectroscopic intensities between the 
state \2b) in the even-even nucleus and all possible states in its odd-even counterpart are 
presented in Table |. In the first column the quantum labels corresponding to the state 
in the odd-even nucleus are specified. Two calculations are shown, the first one under the 
label gs for the traditional operator given in Eq. (||), and the second one under the label 
gs+ for the full operator proposed in this work, Eqs. (^) and (pSj). Important differences 
between both calculations can be seen for some transfers. Let us discuss briefly the transfers 
involving the state \2b) in the even-even nucleus and the states |[2, 0](0, 0)(0, 0)0; 1/2) and 
|[l,l](l,l)(l,l)l;l/2) in the odd-even nucleus. The state |[2, 0](0, 0)(0, 0)0; 1/2) has just 
one component 1(0^ x (see table XIII in Ref. |0), thus the connection with the 

state \2b) in the even-even nucleus can only be done through the terms [s^ da\i^'^^^'^'^/'^'^ in 
the transfer operator. The new terms proposed in this work give no contribution to these 
transfers as it is seen in Table |. The value zero for the j = 5/2 transfer is due to the value 
^5/2 ~ selected. On the other hand, the state |[1, 1](1, 1)(1, 1)1; 1/2) has two components 

(see table XIII in Ref. [H): ^^\{2b x 4/2)^^^^^) + ^[2/^\{2b x 0^/2)^^^^^)- The connection 

with the state \2b) in the even-even nucleus can be done only by the fl]=3/2 5/2 ^"^^ 

terms. The values under the label gs are those coming just from a\. The values under 



the label gs+ include the contributions from aj and those coming from the proposed terms 



[d'^da^). It can be seen that these contributions are important, constructive in one case and 
destructive in the other case. 

It can be in Table | observed that some forbidden transfers with the traditional operator 
given in Eq. (|l]) become allowed with the new one, Eq. (j^). We have performed calculations 
with larger number of bosons with IBA-1 and IBFA-1 and with the two types of bosons 
(neutron and proton) (IBA-2,IBFA-2) and in every case sizeable differences are observed. 

In conclusion, we have calculated extra Aw = 1 terms in the bosonic expansion of the 
transfer operator which connects states of IBA with states of IBFA. Through the comparisons 
between the transfer intensities calculated using the transfer operator with and without these 
terms in the U^^{5) limit of the pseudo-spin symmetry U^^{6)®U^{2) |^T|, we have shown 
that the added terms are not negligible when compared to the ones taken into account up to 
now and expect improvements in the description of the one nucleon transfer intensities using 
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this modified operator. Calculations in the Pt region are in progress and will we published 
somewhere else [p^] . 

The influence of this modified operator for the calculation of beta-decay within the 
interacting boson model remains to be studied. The matrix elements needed in this case are 
of the same kind as the ones needed in order to calculate spectroscopic intensities. However, 
in beta-decay all transfers connecting an initial state in an odd-even nucleus (not necessarily 
the ground state) to all possible states in the even-even intermediate nucleus and then all 
transfers from these intermediate states to the final state (not necessarily the ground state 
too) in the even-odd nucleus have to be calculated. Thus, many transfers contribute and it 
is important to study the changes induced by the new transfer operator in its description. 

Work supported in part by Spanish CICYT under contract PB98-1111. 
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TABLES 



TABLE L Calculated spectroscopic intensities for one-nucleon transfer between an even-even 
nucleus with = 1 and an odd-even U^^{5) nucleus with = 1 plus one fermion. Results 
for transferred j = 1/2, j = 3/2, and j = 5/2 are shown. Calculations made with (|l|) and 
(|5|) are displayed under the labels gs and gs+, respectively. Only transfers involving the state 
\2b) = |[1](1)(1)2) in the even-even core nucleus are included, since those involving the state 
\0b) = I [1](0)(0)0) in the even-even core have no contribution from the added terms to the transfer 
operator. The quantum numbers in the first column specify the state in the odd-even nucleus. 



[Ni,N2]{ni,n2){vi,V2)L;J 
odd-even state 


j 


=1/2 

gs+ 


j 

gs 


=3/2 

gs+ 


j 

gs 


=5/2 

gs+ 


[1,1](1,1)(1,1)1 


1/2 






0.900 


0.477 


0.800 


1.542 


[2,0] (0,0) (0,0)0 


1/2 






0.539 


0.539 


0.000 


0.000 


[2,0] (2,0) (0,0)0 


1/2 






0.600 


0.980 


1.200 


0.844 


[1,1](1,1)(1,1)1 


3/2 


0.000 


0.007 


0.900 


1.638 


2.800 


2.181 


[1,1](1,0)(1,0)2 


3/2 


3.007 


2.526 


0.217 


0.067 


0.000 


0.003 


[2,0](1,0)(1,0)2 


3/2 


0.071 


0.015 


0.217 


0.451 


0.000 


0.003 


[2,0] (2,0) (2, 0)2 


3/2 


0.000 


0.053 


2.100 


1.659 


1.200 


2.164 


[1,1](1,0)(1,0)2 


5/2 


3.546 


5.168 


0.031 


0.220 


0.000 


0.027 


[1,1](1,1)(1,1)3 


5/2 


0.000 


0.142 


3.600 


3.217 


1.200 


2.441 


[2,0](1,0)(1,0)2 


5/2 


0.321 


0.916 


0.031 


0.014 


0.000 


0.027 


[2,0] (2,0) (2, 0)2 


5/2 


0.000 


0.000 


0.900 


1.662 


4.800 


4.921 


[1,1](1,1)(1,1)3 


7/2 






0.600 


0.600 


7.200 


7.200 


[2,0] (2,0) (2, 0)4 


7/2 






5.400 


5.400 


0.800 


0.800 


[2,0] (2,0) (2, 0)4 


9/2 










10.000 


10.000 
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